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Abstract.  Direct  numerical  simulation(DNS)  of  a fully-developed  open- 
channel  flow  has  been  carried  out  for  subcritical  Pi-oude  numbers  at  a fixed 
friction-velocity  Reynolds  number  of  180.  The  free  surface  is  approximated 
by  small-amplitude  wave  theory.  The  results  are  presented  emphasizing  the 
effects  of  the  Froude  number  on  turbulence  quantities  related  to  the  free 
surface  fluctuations.  The  amplitude  of  the  free-surface  fluctuation  increases 
as  the  square  of  the  Froude  number.  Wliile  statistical  quantities  involving 
the  vertical  velocity  component  are  most  influenced  by  the  Froude  number 
but  only  in  the  region  close  to  the  free  surface,  quantities  related  to  the 
pressure  fluctuations  are  influenced  over  much  wider  region. 


1.  Introduction 

Turbulent  flows  in  open-channels  are  important  in  engineering  and  envi- 
ronmental applications,  but  their  understanding  and  prediction  are  not 
as  satisfactory  as  those  for  flows  in  closed  ducts.  Turbulence  mechanisms 
near  a free  surface  have  been  investigated  both  experimentally  and  numer- 
ically(e.g.  Kumar  et  ai,  1998  and  Nagaosa,  1999),  but  no  specific  methods 
of  modeling  in  practical  engineering  calculations  are  established.  It  may  be 
noted  that  the  recent  advancements  of  the  wall  turbulence  modeling  have 
been  made  possible  by  detailed  databases  obtained  by  DNS.  DNS  for  open- 
channel  flows  have  been  reported,  but  most  of  these  simulations  assume 
that  the  free  surface  is  a rigid  slip  surface  and  its  vertical  movement  is  ne- 
glected(e.g.  Nagaosa,  1999).  In  the  DNS  conducted  by  Komori  et  af.(1993), 
the  Pi-oude  number  is  so  small  that  there  are  no  recognizable  differences 
from  the  slip-surface  flow.  Borue  et  af.(1995)  applied  the  linearized  free- 
surface  boundary  conditions  and  obtained  the  dynamics  associated  with 
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the  turbulence/interface  interactions.  They  did  not  provide  the  type  of 
database  that  would  be  required  to  refine  and  validate  turbulence  models 
for  open-channel  flows.  In  the  present  paper,  we  report  the  results  of  DNS 
of  a fully  developed  turbulent  flow  in  a two-dimensional  open  channel  with 
deforming  free  surface  approximated  by  the  small-amplitude  wave  assump- 
tion using  a finite  difference  method.  We  obtain  a detailed  database  of  the 
mean  flow  and  the  turbulence  quantities  for  subcritical  Proude  numbers. 


2.  Basic  Equations  and  Boundary  Conditions 


We  consider  a fully  developed  incompressible  open-channel  flow  over  a flat 
smooth  floor  as  shown  in  Fig.l.  The  governing  equations  are  the  Navier- 
Stokes  equations  and  continuity  equation  for  incompressible  fluid  with  grav- 
itational acceleration  g.  The  boundary  conditions  are  the  no-slip  condition 
on  the  bottom  floor (x2  = 0),  and  the  free-surface  conditions  as  described 
in  the  following.  The  flow  field  is  supposed  to  extend  infinitely  in  the 
strearawise(a;i)  and  the  cross-flow(a;3)  directions  and  the  periodic  bound- 
ary conditions  are  applied  at  the  boundaries  of  the  finite  computational 
domain.  We  denote  the  time  mean  of  instantaneous  quantity  / by  / or  F, 
the  deviation  from  the  mean  by  /,  the  root  mean  square  by  and  the 
value  at  X2  = a by  /la.  The  free  surface  movement  is  assumed  mild  and  it 
can  be  described  by  a smooth  continuous  surface  which  deforms  continu- 
ously so  that  its  position  can  be  represented  by  a single- valued  continuous 
function  by  X2  — h{xi,  x^,  t).  Function  h satisfies  the  kinematic  condition 


dh 

dt 


dh 


dh 


(1) 


where  Uj  is  the  instantaneous  velocity  .component  in  Xj  direction.  If  we 
assume  that  the  fluid  above  the  free  surface  has  negligible  density  and  the 
pressure  is  constant  zero,  and  that  the  surface  tension  can  also  be  neglected, 
the  dynamic  conditions  are  that  both  normal  and  tangential  stresses  on  the 
instantaneous  free  surface  are  zero.  If  the  displacement  of  the  instantaneous 
free  surface  from  the  mean  position  and  the  slopes  of  the  free  surface  is 
small,  the  quantities  on  it  may  be  replaced  by  the  values  on  the  mean 
position  and  the  normal  and  tangential  directions  may  be  replaced  by  the 
vertical  and  horizontal  directions,  respectively,  then  the  first  approximation 
gives 
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where  P = p — g2h,  p is  the  instantaneous  pressure  divided  by  the  density, 
Q2  is  the  component  of  gravitational  acceleration  in  X2  direction  and  u is 
the  kinematic  viscosity.  Since  these  boundary  conditions  are  all  in  terms  of 
the  quantities  at  the  mean  position  of  the  free  surface,  the  flow  calculations 
can  be  made  with  the  fixed  grid  in  the  region  X2  < H. 


3.  Numerical  Methods 

The  governing  equations  are  solved  by  a finite  difference  technique  based 
on  the  SMAC  method  on  a Cartesian  staggered  grid.  The  spatial  deriva- 
tives are  descretized  using  the  second-order  conservative  difference  scheme 
and  time  advancing  is  done  by  the  second-order  Adams-Bashforth  method. 
Tlie  position  of  the  free  surface  h is  solved  by  descretizing  the  spatial 
derivative  terms  in  Eq.(2)  using  the  fifth-order  upwind-shifted  interpola- 
tioii(USI)  scheme(Kajishima,  1994)  and  time  advancing  by  the  third-order 
Adams-Bashforth  method.  The  Poisson  equation  for  pressure  is  solved  by 
Fourier-transforming  the  equations  in  the  horizontal  directions  and  the  re- 
sulting equation  is  solved  by  the  tri-diagonal  algorithm.  The  calculations 
were  conducted  for  the  Reynolds  number  Rer  based  on  the  friction  ve- 
locity Ut  and  the  mean  flow  depth  H of  180,  which  is  the  same  as  the 
closed-channel  flow  simulation  of  Moser  et  a/.(1999)(MKM)  and  Yokojima 
& Nakayama(2000)(YN)  for  a similar  flow  with  free-slip  boundary  condition 
(hereafter  referred  to  as  the  ’slip  channel’).  The  Proude  numbers  based  on 
the  average  velocity  and  the  flow  depth  are  0.3,  0.6  and  0.9.  The  number  of 
grid  points  is  128^;  the  spatial  resolutions  are  = 9,  AxJ  = 0.27 - 2.71, 
A^3  = 4.5.  Here  the  superscript  refers  to  the  nondimensionalized  quan- 
tities by  Ur  and  i/.  With  a time  increment  of  At'^  = 0.0114,  about  200,000 
time  steps  are  used  to  obtain  the  well-converged  turbulence  quantities  av- 
eraged over  space  and  time. 


4.  Results  and  Discussion 

Fig.2  shows  properties  of  the  calculated  instantaneous  flow  field  for  the  slip- 
channel  flow  at  Fr  = 0 and  the  present  fluctuating  free  surface  computation 
for  Fr  = 0.6.  The  figures  on  the  top  are  the  samples  of  instantaneous  flow 
represented  by  the  iso-surfaces  of  the  streamwise  velocity  u\.  The  instan- 
taneous position  of  the  free  surface  h and  the  velocity  vectors  projected  in 
the  mean  free-surface  position  are  shown  in  the  middle  and  at  the  bottom, 
respectively.  The  well-known  streak  structure  of  the  wall  turbulence  is  seen 
in  both  cases  near  the  bottom  floor  and  the  overall  simulation  method  ap- 
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pears  satisfactory.  In  Fr  = 0 case,  h is  not  available  and  the  instantaneous 
pressure  at  the  free  surface  P\h  is  shown  instead.  The  surface  pressure  and 
h appear  to  show  very  a similar  trend.  The  free  surface  fluctuation  and  the 
velocity  vectors  do  not  seem  to  have  any  prefered  directionality  and  their 
scales  appear  to  be  much  larger  than  the  viscous  scales  of  the  near-wall  flow. 
These  tendencies  agree  with  the  results  of  Thomas  k Williams(1995)(TW) 
who  reported  preliminary  results  of  a direct  simulation  using  ’volume  of 
fluid’  method  to  track  the  movement  of  the  free  surface.  The  velocity  vec- 
tors on  the  free  surface  indicate  that  the  flow  is  generally  divergent  in  the 
horizontal  plane  where  the  upwelling  bulges  with  positive  h are  seen  and 
convergent  in  the  area  the  free  surface  is  lower  than  the  average  position. 

The  profiles  of  computed  RMS  surface  fluctuation  are  compared  with 
experimental  results  of  Nakayama  et  al.(2000)(EXPN00)  and  Nakayama 
(1997)(EXPN97)  in  Fig. 3.  The  Reynolds  numbers  of  the  experiments  are 
indicated  in  terms  of  the  bulk  Reynolds  number  defined  by  the  average 
velocity  and  the  mean  depth  H.  It  shows  that  the  present  DNS  results 
are  a little  lower  than  the  measurements  but  show  the  same  increasing 
trends  with  respect  to  Fr  and  are  very  close  to  the  results  of  TW.  Though 
not  shown,  the  computed  skewness  and  flatness  factors  of  h also  agree 
well  with  those  of  measurements,  which  indicates  that  the  small-amplitude 
approximation  used  here  is  appropriate  for  the  DNS  of  subcritical  flows. 

Fig.4  shows  the  computed  mean  velocity  profiles  compared  with  the 
closed-channel  and  slip-channel  flow  simulations.  It  shows  that  the  velocity 
profile  of  the  open-channel  flow  follows  the  logarithmic  distribution  up  to 
a point  very  close  to  the  free  surface  and  it  is  not  influenced  very  much 
by  the  free-surface  movements.  The  distributions  of  the  RMS  fluctuations 
of  velocity  components  and  the  Reynolds  shear  stress  are  shown  in  Fig.5. 
It  is  seen  that  the  distributions  of  all  stress  components  are  very  close  to 
the  closed-channel  results  of  Moser  et  al.(1999)  over  most  of  the  channel. 
The  open-channel  flow  results  deviate  from  the  closed  channel  results  in 
the  region  120  < xj  < 180.  Fig.6  is  an  enlarged  plot  of  n™®  which  is 
the  quantity  most  sensitively  affected  by  the  free  surface  and  Fr. 
increases  with  roughly  square  of  Fr,  which  is  the  same  as  Fig. 7 is  a 

plot  of  the  RMS  pressure  fluctuation.  Unlike  the  velocity  fluctuations, 
is  seen  to  be  influenced  by  the  free  surface  and  the  value  of  Fr  throughout 
the  entire  channel  depth.  Since  the  Reynolds  stress  component  is  most 
influenced  by  the  Froude  number,  the  terms  in  the  transport  equation  for 
this  stress  are  shown  in  Fig.8.  Fig.8(a)  is  the  distribution  across  the  entire 
channel  and  Figs.8(b)-(d)  show  enlarged  plots  near  the  free  surface  for 
Fr  = 0.3,  0.6  and  0.9,  respectively,  each  compared  with  Fr  ~ 0 case.  It  is 
seen  that  near  ^e  free  surface  the  pressure-strain  term,  which  works  like  a 
generation  for  drops  and  the  pressure  diffusion  increases.  It  is  further 
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seen  in  the  enlarged  plots  of  Figs.8(b)-(d)  that  the  pressure  diffusion  term 
changes  sharply  from  loss  to  gain  within  about  20  viscous  units  from  the 
free  surface.  Its  distribution  and  that  of  the  turbulent  diffusion  appear  to 
be  influenced  by  Fr  over  at  least  upper  half  of  the  channel. 

5.  Conclusions 

Direct  numerical  simulation  of  a fully  developed  turbulent  flow  in  a two- 
dimensional  open  channel  has  been  successfully  performed.  The  present 
results  with  the  assumption  of  small  displacement  of  the  free  surface  from 
its  mean  position  appear  to  represent  the  open-channel  flow  with  subcritical 
Pi'oude  numbers  properly,  and  show  the  effects  due  to  the  moving  free 
surface.  Most  of  the  quantities  including  the  mean  velocity  profiles  and 
the  Reynolds  stresses  in  the  horizontal  directions  are  not  influenced  by  the 
free-surface  movements.  Tlie  vertical  velocity  fluctuation  and  the  pressure 
fluctuation  are  the  quantities  influenced  by  the  free  surface.  The  vertical 
fluctuation  near  the  free  surface  increases  with  Fr  in  the  region  very  close 
to  the  free  surface,  while  the  pressure  fluctuation  is  influenced  over  almost 
entire  depth  of  the  channel.  This  influence  appears  to  be  effected  by  the 
readjustment  by  the  pressure  diffusion. 
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(a)  Low-speed  streaks,  Fr  = 0 


(b)  Low-speed  streaks,  Fr  = 0.6 


(e)  Projected  velocity  vectors,  Fr  = 0 (f)  Projected  velocity  vectors,  Fr  = 0.6 

Figure  2.  Instantaneous  flow  features. 


